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The paper discusses the transformation of decorated Ising models into an effective undecorated 
spin models, using the most general Hamiltonian for interacting Ising models including a long range 
and high order interactions. The inverse of a Vandermonde matrix with equidistant nodes [— s,s] 
is used to obtain an analytical expression of the transformation. This kind of transformation is 
very useful to obtain the partition function of decorated systems. The method presented by Fisher 
is also extended, in order to obtain the correlation functions of the decorated Ising models trans- 
forming into an effective undecorated Ising models. We apply this transformation to a particular 
mixed spin-(l/2,l) and (1/2,2) square lattice with only nearest site interaction. This model could 
be transformed into an effective uniform spin-S' square lattice with nearest and next-nearest interac- 
tion, furthermore the effective Hamiltonian also include combinations of three-body and four-body 
interactions, particularly we considered spin 1 and 2. 

1. INTRODUCTION 

Two-leg transformation was first introduced by Fisher[l| in the 1950's decade, after which this transformation was 
used widely for one-dimensional decorated models, such as discussed recently for the tetramer Ising-Heisenberg bond- 
alternating chain as a model^stem for Cu(3-Chloropyridine)2(N3)2 in reference Another decorated model has 
also been considered recently[3|, one that can be applied even for two-dimensional decorated Ising modelsQ or higher 
dimensions. 

The Ising models with multi-spin interactions have been extensively investigated both theoretically and experimen- 
tally. In this sense, we present a review of the start-triangle transformation[l|. This kind of transformation also was 
discussed in detail by SyoziQ. This transformation is widely used for two-dimensional models, the most well known of 
these being the Ising kagome lattice[^, the honeycomb latticeQ, and other two-dimensional models like Ising model 
on pentagonal lattice[l|. 

Another illustration of the application of this transformation was performed for the Ising-Heisenberg diamond 
chain[3,i|. This transformation can also be applied to higher dimension decorated lattice such as, two-dimensional 
decorated Ising-Heisenberg modelsQ and two-dimensional doubly decorated Ising-Heisenberg models(lo|. 

Recently several real systems motivate to investigate in this kinds of transformation, such as the discovered two- 
dimensional magnetic materials CU9X2 (cpa)e ■XH2O (cpa=2-carboxypentonic acid; X=F, CI, Br) where the Cu spins 
stands on the triangular kagome lattice[lj| with Heisenberg interaction type. Liquid crystals networks composed by 
pentagonal, square and triangular cylinders Other recent investigation about the crystal structure of solvated 
[Zn(tpt)2/3(SiF,)(H20)2- (MeOH)] [tpt = 2,4,6-tris(4-pyridyl)-l,3,5-triazine] networks with the (10,3)-a topology[l3|. 

This paper is organized as follows. In section 2, we start with a review of star-triangle transformation[l| where we 
comment their possible extension, afterwards we extend it formally to the general case of the m-leg spin-1/2 system. 
In section 3 we present the transformation for the general case of the m-leg spin-S", using matrix formulation. In 
Section 4 we discuss the correlation function for the general case of the m-leg spin-S" system also using the matrix 
formulation. In section 5, this transformation is applied to the mixed spin-(l/2,S') of square Ising lattice with only 
nearest interaction, spite this model cannot be mapped into exactly solvable model, we even could discuss their 
critical point behavior. Using a double transformation we study the spin-S* square lattice with nearest and next- 
nearest interaction, furthermore the effective Hamiltonian also include combinations of three-body and four-body 
interactions. Finally in section 6 we present our conclusion. 

2. TRANSFORMATION OF DECORATED SPIN-1/2 MODELS 

2.1. Star-triangle transformation 

The star- Ising model with an arbitrary decorated spin Sq, and the spin-1/2 of legs cri, (T2 and 0-3 is presented in 
the figure 1. Conveniently the Hamiltonian will be defined from now in units oi (3 — 1/kT, being k the Boltzmann 
constant and T the absolute temperature. Therefore the Hamiltonian is given by 
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+ •^i^o,i'^"'^i'^3 + Jo^liSo<y20-3 + j[^liSacria2(J3, (1) 
here the parameter JqI) o represents the magnetic field acting on spin Sq, J['^qq, Jq^I q and Jq^q i are the exchange 

f3) (3) 

parameter between Sq and cti, 0-2 and (T3 respectively. The parameters of the three-spin interaction are Jj; / g, 5 and 

(3) 

Jg 1 1 which corresponds to the three spins products So<7ia2, Soaia^ and S'oCT20'3 respectively. Finally the interaction 

(3) 

parameter Jl ii corresponds to the four-body interaction S'ocricr20'3- It is important to notice that the first fourth 
terms of the Hamiltonian already was considered early by Fisher 

Now the Hamiltonian of the triangle- Ising model represented in Fig[T] (right) can also be expressed, in a very general 
way, as follow 

(72, CF3) ==Jofi,0 + Jhlo'^i + Jo!lo'^2 + Jo^liCTs + j[^lo'7l(^2 + 

I 7(3) I 7(3) , 7(3) /„N 

+ >^1,0,1^1'''3 + >^04,1^2Cr3 + ^1 1 ^CTiCT2Cr3, [Zj 

— (3) 

where Jq g g correspond to the effective parameter corresponding to a constant energy, the coefficients of a; (with 
i = 1..3) represents the effective parameter of a magnetic field, the coefficients of aiaj (with = 1..3) are 

the standard bilinear coupHng effective parameter, and the coefficient of (Ti(T20'3 represent the effective three-linear 
interaction. 

Besides adding a magnetic field parameter to the standard star-triangle transformation proposed by Fisher[l|, we 
also included the three-linear parameter, in order to solve the algebraic system equation consistently. We can verify 
that there are eight equations and eight parameters to be obtained. 

Carrying out a partial trace over the variable of the decorating system, it is reduced to an effective partition function 
with spin-1/2 Ising model. Therefore we write down the partition function for a given decorated system 

Z(/3) = e^''^''«.«Z(/3), (3) 

where Nd is the number of decorations of the lattice. Considering a partial summation on decorated particles (Sq), 
we have the following amount, which we will call the associated Boltzmann weight 

W'(^n'Ti,^2,a3) =tr5„(e«'''(^-'^---'^^)). (4) 

On the other hand the Boltzmann factor M^'-^-'((7i, 0-2, 173) in the transformed (undecorated) Hamiltonian, becomes 
W<'^\ai,a2,a3) = exp ( J^'g^ g + j[%^ai + J^%^C72 + J^J^a^ + ji%o<Ji<72 + 

+ -'lio,l^l'^3 + l'^2(^3 + •^a.lf^l^2Cr3). (5) 

The transformation of the parameter of Jni.712 ,n3 into a function of Jni,n2 ,ti3 is obtained, relating the W^^'^ (cri , 0-2 , cr3) = 
, o'2,o'3)- Therefore from eqs.lH]) and ([5]), we get the solution of the new parameters as a function of the 
decorated Hamiltonian parameters, given by 

J^%.,ns = l E <'^r<^r^^iW^'H<^U'^2,'J3)), (6) 
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where = 0, 1 with i = 1,2, 3. 

K we consider the arbitrary decorated spin Sq as a simple spin-1/2, this transformation is known as duahty of the 
star-triangle relation^ and is useful to solve other spin models such as kagome like models. 

The solution showed in eq.|[6]) corresponds to the eight solutions. A particular case of this solution becomes the 

solution obtained by Fisher[l|, when we consider Jpg q = Jq\' q = Jq^q ^ = and i = 0, leaving thus only four free 
parameters to be determined instead of eight. 

This particular solution is widely applied to a large variety of two dimensional lattice. In a recent paper, Lackova 
et aZ.[lJ| discussed an exact results of a mixed spin-1/2 and spin-1 transverse Ising model with two- and four-spin 
interactions and crystal field on the honeycomb lattice. 



2.2. m-leg star-polygon transformation 

The transformation of the decorated model presented by Fisher[l| could also be extended to the cross-square 
transformation as considered recently by Savvidy[ll|, to study the systems with exponentially degenerated vacuum 
state. A particular case of this transformation was discussed in reference with the bilinear and the four-linear 
terms as the non-null parameters, thus maintaining only three free parameters in the Hamiltonian, which is necessary 
to study the gonihedric model|16l|. 

Motivated by the previous result, we can extended the formulation to a general case, where we consider m particles 
interacting with a central spin Sq. This kind of model can be used to study two or three dimensional lattices with 
high order coordination number. The Hamiltonian of the decorated model with non-linear interaction and central 
spin Sq can be expressed as 




FIG. 2: Schematic representation of m-leg decorated Ising model (left) and m-side polygon Ising model (right) 



{ni}„=04 i=l 

where >/|™}^ is the parameter interaction of 01^ i ''"f' '^ith sub-indexes {ai}m- We denote the set {<7i}m = 
{(Ji, . . . , dm}, similarly by {ni}m = {ni, . . . , n,„}, whereas the super indexes are related to the number of legs or 
coordination number of decorated spin. 

On the other hand, the undecorated Hamiltonian into which the decorated Hamiltonian ((7| will be transformed, is 
given by 

m 

H("'({aj„o= E ^{:in-r- (8) 

{ni}„,=0,l i=l 

The transformation of the parameters J^n\ ^ function of J^n\ ^ire obtained in a similar way as in the previous 
case. Writing the associated Boltzmann weight, as a trace over the spin-iSo of the decorated system, 

M/(™)({^^}„^) = tr5o(c«'-''(^-{->")). (9) 
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On the other hand the Boltzmann weight function W^"^\{ai}m.), for the transformed Hamiltonian can be expressed 

as 

m 

w("'\w^u)= n ^-p(^{:m-ro- (10) 

{cri}m = ±l «=1 

Assummff that - W^"'H{<J,}^) and eq.© are known, we substitute them in eq. lfTO|) . getting the solution 

of the unknown parameters which reads as 

4:1 = ^ E iUO^-i^'-Hw^u)). (11) 

{<T,}„=±1 i=l 

Rewritten the eq. lfTTj) a httle different, we have 

jt^L = E .9{n,}„({aJ™)ln(M^M({a,}„0), (12) 

where the function g{„;}^ depends only of spins ai, 

^ m 

9{n,U ^ n "^'"^ ' (1^) 

i=l 

with Ui = 0,1 and i = 1, 2, . . . , to. Eq. (fT2l) corresponds to 2™ solutions of the unknown parameters • 

Let us consider a special case as an example, without losing its general properties, suppose we consider the central 
spin So = ±1/2, then the eq.© is reduced to 

m 

WM({aO™) = 2cosh( J{nLllO- (14) 

3. TRANSFORMATION FOR SPIN-s MODELS 

One interesting extension considered here is the higher spins of decorated models, then for this purpose we first 
consider the two-leg spin-1 model, where we introduce a matrix formaHsm to simplify our notation. After that we 
will discuss the fully general spin-s m-leg star-polygon transformation. 

3.1. Two-leg transformation for spin-1 model 

In this section we discuss the spin-1 case. Recently similar situation was discussed by Fireman et a/.(l8l| where 
they consider the Ising model with bi-Hnear, bi-quadratic, single ion anisotropy and Zeeman interaction, which was 
mapped into an effective Blume Emery Griffiths (BEG) model[l^. Extending to the spin-1 model, it is possible to 
write the Hamiltonian according to the definition in eq.©, 

2 

n^'HSo,s,,s2)= ^ 4%,SosrsT, (15) 

ni ,712—0 

(2) 

where Jni,n2 corresponds to interaction parameter, here si and S2 represents the spin-1, whereas S'o could be any 
mechanical system. The effective Hamiltonian of the undecorated system also can be expressed in analogy with eq.®, 
which would give us, 

■H^^\s,,s,)= JL%24'4', (16) 

Til ,712—0 



5 



with JAi.ria being the interaction parameters of the transformed system. Considering a partial summation of the 
partition function where involves the summation (trace) on spin Sq, we obtain as before the so called associated 
Boltzmann weight, given by, 



(17) 



To solve the unknown parameter we introduce the matrix formalism which is an appropriate representation to express 
these equations in terms of the known parameters. Using the matrix notation we write the eq. (fT5| . simply as follow 



H(2)(5o) = 5oV(1'1)j(2), 
where V^i-i) and J^^) are defined respectively by 



(18) 



and J(2) 



(19) 
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To transform the decorated Hamiltonian in eq. ljlSp into an undecorated Hamiltonian, as given in eq. lfT6|) we find 
the parameters of effective Hamiltonian as a function of the original decorated Hamiltonian. Therefore, we define the 
associated Boltzmann weight function W'^^^si, S2) as 
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(20) 



To make a complete matrix representation, we define the following function 7?.'^'(si, S2) = ln(M^(^)(si, S2)). The 
function T?.'^' (si, S2) can be expressed as a vector R^^^. Using the matrix notation and considering VF*-^^ (si, S2) = 
W^^\si, S2), we have the following equation 



r(2) ^ v(^'^)j(2 



(21) 



The solution of the equation above could be obtained by taking the inverse of the matrix V'-^'^-', denoted from now 
on by V^^'^' — (V'^'^^)^^, The solution is then expressed by 



j(2) ^ y(1.1)r(2)_ 



(22) 



But the V^^'^) can be written also as V^^'-^) = V^^^ (8) V^^) and its inverse becomes V^^'^^ = V^^^ V^^^. In the present 
case the matrix V^^^ and V^^^ could be written expHcitly as 




and V(i) 




(23) 



Certainly it is simpler to evaluate the inverse of the reduced matrix V^^^ instead of the large matrix V'^^^^. 

As an example we can consider the non-uniform two-leg Ising spin transformation with central spin Sq ~ 1, without 
losing any generality. For this particular case the eq. (flT)) will be expressed as 



M^('Hsi,^2)=l + 2cosh( 

7i\ ,n2— 



(24) 



We remark that this transformation is carried out for every unitary cell and can be appHed to the one-dimensional 
01, two-dimensional [l3| and other high dimensional spin lattices. 
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3.2. Two-leg transformation for spin-s model 

For a higher spin the Hamiltonian can be written in a general way using the previous matrix notation 

H(2)(5'o) = S'oV('''^'j(2) = S-qV^") ® V(^'J(2), (25) 

where the dimension of the matrix V^'''*^ is (2s + 1)^ x (2s + 1)^, meanwhile, the dimension of the matrix V'^'*^ is 
(2s + 1) X (2s + 1). We can see that the matrix V^'') is a Vandermonde matrix with equidistant nodes [— s, s], and the 
elements of the nodes are xj which corresponds only to the magnetic moments of the spin-s. The elements can be 
appropriately expressed as Xj = — s — 1 + j, with j = 1, 2, . . . , 2s + 1. The explicit representation of the Vandermonde 
matrix is given by 



xi 

X2 



1 X2s 
\l X2S+1 X^ 



"-23 



Xi 



X2s 



2s+l -^23+1 



4' \ 

X2 



T.2S 

^2s 



(26) 



Using he matrix notation, the function R(^^(si,S2) is defined in a similar way as was defined in eq.((2T|), therefore 
we have a vector function given by 



R(2)=ln(W(2)) = ln(trs„(cH<''(S"))). 



(27) 



The effective Hamiltonian of the undecorated system also can be written in analogy to eq. (|T6l) . Thus we have 

H(2) = v(-^'^)j(2) = V(^) ® v(^)j(2). (28) 

The Hamiltonian considered in eq. l(28|) . can contain at most (2s + 1)^ parameters. This also means that the 
dimension of the vectors J^^' and R^^' is (2s + 1)^. Taking the inverse of the matrix V^''''*^ we are able to express the 
new parameters as a function of the known parameters defined by eq. l(25|) , the new parameters then are the elements 
of 



j(2) ^ y(s,s) p(2) ^ y{s) ^ y{s) p,{2) ^ 



(29) 



The inverse of the matrix V''*) could be solved using the recursive equation presented recently by Eisinberg et a/.(20l|. 
where is discussed an generic algorithm to obtain the elements of the inverse of the Vandermonde matrix V^''). 
Therefore following that algorithmfl^l, we find that the elements of the matrix V''') could be written explicitly as 



V 



is) (-!)'+■' 

(2s + l-j)!0--l) 



1) 



k — i 



2s + 2 
k + 1 



^t+l \^ s+l) 



(30) 



where 



represents the first kind of Stirling number, whereas F,^^^ represents the hyper-geometric function as 



defined in 

It could be interesting to normalize the spin-s. In these cases, it is possible to rewrite the elements of the Vander- 
monde matrix V^") with equidistant nodes [—1, 1] , whose elements are given by Sj — (— s — 1 4- j)/s. The solution of 
this case has already be found by Eisinberg et a/.[20|- Here we present the same solution but using the hyper-geometric 
function explicitly, as follows 



2s + l 



^■^ (2s + l-j)!(j-l) 



fe=l 



2s-f 2 
k + 1 



s+l 



(31) 



Certainly this solution could be useful when the spin-s is large. 



7 



3.3. The m-leg star-polygon transformation for arbitrary spin 

This section presents the general extension for the m-leg star-polygon transformation of decorated Ising spin model. 
Based on the previous results, the considered Hamiltonian will be written as follows 

H(™H^o) = (g) V^^'^^oJ''"^ (32) 

1=1 

where Si corresponds to an arbitrary spin for each leg. Assuming that we have m legs of the star or m edge of the 
polygon with 5*0 being the central spin, then the dimension of the matrix is (2si + l) x (2si + l). The Hamiltonian 
that we consider can contain at most Ililil^Si + 1) parameters, which means that the dimension of the vector J^™' 
and R^™^ is n"=i(2si + !)• The elements of the vector J^™) are the parameters of the original system. 
We define in analogy to the eq. (|27)) . a vector function for the following amount 

R(")=ln(tr5o(e^'"'(^°0). (33) 
The effective Hamiltonian in an undecorated system also can be expressed in analogy to (f28ll . thus we have 

m 

H(") =(g)V^''')j(™), (34) 

■i=i 

where the J^™^ is a vector whose elements are the parameters of the transformed system. Therefore, the new parameters 
could be written using the vector representation as 

m 

J(™) = (g)V(''')R(™), (35) 

i=l 

we remark that the matrix V*^**'^ is the inverse of the matrix V'^'**). The corresponding Boltzmann weigth could be 
written as follow 



M/({s„})=e««^-». (36) 

The explicit form of matrix V''*' and its respective inverse is presented in appendix A, up to spin-3. For an arbitrary 
spin-s, the matrix V^**^ is given by eq. l(26|) and the elements of the matrix inverse are given by the eq. (|30| . 



4. THE CORRELATION FUNCTION OF DECORATED SPIN MODELS 



The correlation function in the transformed Hamiltonian involving the central spin 5*0 can be obtained in a very 
similar way such as initially proposed by Fisher[l|. Here we present a very general extension to obtain the correlation 
function, using the correlation function of effective systems (undecorated systems). 

The correlation function involving spins such as So, Ski, ■ ■ ■ , Sk^, with Sk^ any arbitrary spin of system, can be 
written as 

{SoSk,...Sk^) = ltrs„,{,,j(5oSfe,...Sfc^e«""'(^-^--^")), (37) 

where the trace is performed over all 5*0 and {si} — {si,S2, ■ ■ ■ ,Sm}, whereas Z is the partition function for the m-leg 
Ising model. For simplicity, we consider, as before the Hamiltonian in units oi (3 = — 1/fcT. The partial summation 
over So can be expressed as follows 

C(5i,...,s„0 = E^o^''""'^^"'"-''"^- (38) 

So 

On the other hand we want to write the correlation function (|37| as a function of all spins of the m-leg lattice; for 
that purpose, we can try to represent the C(si, . . . , s,„) as a linear combination of si, . . . , Sm, which reads as 

m 

C{si,...,Sm) = E a(ni,...,n„)]^s"', (39) 
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where Q;(ni, . . . , n™) are coefficients to be determined. Matrix notation is usually very convenient, as was shown 
in the previous section. Then we write the eq. (|39| using matrix notation, which reduces it simply to the following 
expression 

m 

C^'"' = 0V('*')a("), (40) 

i=l 

where the dimension of vectors C^™^ and a^™^ is n"=i(2si + 1)- The elements of the vector C^™^ are given by eq. ((38l) . 
whereas the elements of vector a^™) are the unknown coefficients a(ni, . . . , n„i) to be determined in a similar way as 
was determined by eq. ((35l) . Therefore the solution of the unknown elements of the vector a^™^ is given by 

m 

a(™) =(g)V('^')C(™). (41) 

1=1 

To obtain the explicit form of the elements of a^™) , the inverse of the same Vandermonde matrix V^'') must be 
evaluated for each coordination number. Thus, the correlation function can be written as a linear combination of the 
correlation function of the transformed Hamiltonian, and could be written explicitly as 

m 

(5oSfei . . . Sfe,,) = ^ Q;(ni,...,n„)(.Sfej ...Sfc,. ]^s"'). (42) 

Finally, we are able to write the correlation function including the decorated spin, which can always be expressed 
just as the correlation function of the effective undecorated Ising system. 



5. CRITICAL POINTS OF 2D ISING SPIN-S MODEL WITH UP TO QUARTIC INTERACTION 



The two dimensional lattice is one of the most interesting subject in statistical physics, both experimentallv|l{ 
and theoretically. Several approximation methods, such as the mean-field-t heor v 1 1 9l. 1 23l| . the Bethe approximation|24 
the correlated effective field theoryflBl, the renormalization groupjH], the series expansion methods [27|. the Monte 
Carlo methods [2^1 and the cluster variation methods are used to investigate this interesting lattice. However, an exact 
solution has been obtained only in a very limited cases, mainly in the honeycomb lattices |29l[30| . Some exact results 
with restricted parameters has been investigated by Mi and YangfsH using a non-one-to-one transformation[30|. 

As an example we apply the transformation method presented in the previous section to a lattice with mixed 
spin-1/2 and spin-S. The schematic transformation is displayed in figIS] 



5.1. The mixed Ising spin-(l/2,S) lattice 



It is possible to transform a mixed spin lattice into an effective spin-1/2 lattice, such as presented in the literature 
[s^ l- If the spin-S" site is considered as a decoration of the lattice £, then the Hamiltonian is given by 

= E ^'^^^^ + E (43) 

where K being the first neighbor interaction and D the single ion-anisotropy parameter. The case when D — has 
already been discussed by Tang[33|. With ai we represent the spin-1/2 particle, while with Si we represent the spin-S 
particle. 

We have the following Hamiltonian for the effective spin-1/2 Ising model, which is given by 



'^1/2 = + E "^^O'lCTj + E Ji'^n'^i-i'^ii'^ti^ (44) 

all square 

where Jo means the constant energy, J2 corresponds to the interaction parameter with its summation (i, j) running 
over all the pairs of spins interactions for each unitary cell of the lattice and J4 represents the interaction parameter 
among all spins for a given unitary cell of squared lattice considered. 
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D 



FIG. 3: Schematic representation of square-type lattice {C) with spin-1/2, mixed spin-(l/2,S) on square lattice. 



could be obtain using the equation (|36|) as follow 
Due to the symmetry of the lattice, the Boltzmann weights of the mixed spin-(l/2,S') lattice are related by 



The Boltzmann weight of the eight- vertex model [ 

Wi ^^(2; 2' 2' 2)' ^2 — ^{21 ~i 2' ~)' 



W2 =W3 

W5 =we 



W4, 



(45) 



(46) 
(47) 



where 



Wi 



l + 2E„=iCOsh(2nii:)e 



2En=i cosh(n/i)e^ 



S = integral 
S = half-odd-integral 



and 



W2 



l + 2Eti 



S+ 



2 „«^r>. 



S = integral 
S = half-odd-integral 



W5 



1 + 2 ELi cosh(n/C)e^ 



2EScosh(fif)e»^«; 



5* = integral 
S = half-odd-integral 



(49) 



(50) 



From eqs. (|48|) . (|49|) and ((50|) we can verify the following relation for the Boltzmann factors, wi > > W2, 
for any arbitrary values of K and D. The transformation of the mixed Hamiltonian (|43|) into the effective spin-1/2 
Hamiltonian l(44| . relate their parameters, using the Boltzmann weight function, which reads as 



^0 =g In (wlwjwi) , Jo > 0, 



1 



8 \W2 

8 V 



J2 > 0, 

, J4 e 



(51) 
(52) 

(53) 
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The eight- vertex model has been solved approximately by Fan and Wu [s^l when A/w^-^ <^ 1 with Wm = 
max(u>i, u'2, wa, u'4), and with 

A = ujii02 + — uj^ujq — ojjOJs- (54) 

Unfortunately the free-fermion condition (A = 0), does not satisfy our transformed lattice, unless when D zLco, or 
T — > cx). But even so, we can discuss their critical points behavior, using the critical condition 

wi = W2 + W3 + W4, where W2 = W2 — A/wi. (55) 

We have calculated the critical coupling Dc as a function of the parameter Kc, the curve is displayed in fig. 4, for 
spin-(l/2,l) and spin-(l/2,2). Using the equation ((SS)) in all regions of our plots, we verified the amount jAj/uiJ < 1, 
guarantying the convergence of the approximation. A similar situation was considered by Tang [33 |. 
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FIG. 4: We display the critical regions for Dc as a function of Kc in the case of the free-fermion aproximation. The solid line 
corresponds to spin-1, while with the dashed line we represent spin-2 lattice. 

We can verify that the limiting case of the critical function: when Kc ±00 we have a linear behavior of the 
parameter Dc = —2\Kc\ — ln(2)/2, whereas when Dc — » 00, the parameter Kc ±ln(l + \/2)/S. 

Another possibility of the exact solvable condition could be the symmetric vertex model, but unfortunately this 
model also does not satisfy the symmetric vertex model condition because we have always wi > W2 instead wi =1112. 

We point out that the matrix method calculation of the parameters becomes relevant for higer spin values, the 
method used in this case simplify effectively very tedious calulations. For example, for spin-1 the calculation of the 
interaction parameters, without using the matrix method, even could be obtained using the transformation of reference 
[l|. However for spin-2 case it turns on a very cumbersome calculation, we can obtain easily these parameters using 
the matrix method and algebraic software. 



11 



6. CONCLUSIONS 

This paper presents an extension of the transformation of decorated spin models presented by Fisher[l| into an 
effective Ising-type models, including high order and long range interaction terms for an arbitrary number of particles 
and with an arbitrary spin-s. We have shown that the general transformation of decorated spin models is reduced 
to the calculation of the inverse of the especial case of the Vandermonde matrix with equidistant nodes [— s, s]. The 
matrix formulation of decorated spin models transformation could be manipulated more easily when is studied an 
involving decorated systems. This matrix formulation is the main achievement of this work. 

The correlation function of decorated spin models could also be transformed into an effective Ising-type model, 
using the matrix formulation considered above. We verify that, to obtain the correlation function, it is necessary to 
perform the inverse of the same Vandermonde matrix with equidistant nodes. 

As an application, we present a square-type Ising lattice with mixed spin-1/2 and spin-S", following we transform 
this model into an effective Ising model with spin-1/2. Similar models can involve more complicated calculations in 
the case of higer spin values where the matrix method proposed in our manuscript turns relevant. In this sense we 
conclude that the matrix procedure for obtaining the unknown interaction parameters simplify effectively very tedious 
calculations. Finally, as illustration we give in the appendix expHcitly the Vandermonde matrix up to spin-3. 
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APPENDIX A: SOME VALUES OF THE MATRIX V^^' AND V<^' 



In this appendix we present the matrix V^'') and its inverse V^'*-', for spins up to spin-3. For spin-3/2 
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/I 
1 
1 
1 
1 



5 


25 


125 625 


3125 


I 


i 






! 


f 


f ¥ 




2 


! 


i« ¥ 






i 








I 


1§5 6% 


3f§5 




4 


8 16 


32 



25 75 75 25 3 \ 

25fi 128 128 256 25fi 



3 


25 


75 


25g 




1% 


~6|0 


m 


ff 






if 


f 


f 


¥ 


48^ 


i^^ 




120 


24 


12 



2^6 
9^6 



ii 

-^4 



¥ 



120 



(A3) 



For spin-3 
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